Generated by Cython 3.0.2
Yellow lines hint at Python interaction.
Click on a line that starts with a "+" to see the C code that Cython generated for it.
Raw output: edges.c
+01: # cython: infer_types=True, wraparound=False, nonecheck=False, boundscheck=False, cdivision=True, language_level=3, profile=True, autogen_pxd=True
__pyx_t_4 = __Pyx_PyDict_NewPresized(0); if (unlikely(!__pyx_t_4)) __PYX_ERR(0, 1, __pyx_L1_error) __Pyx_GOTREF(__pyx_t_4); if (PyDict_SetItem(__pyx_d, __pyx_n_s_test, __pyx_t_4) < 0) __PYX_ERR(0, 1, __pyx_L1_error) __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
02:
03: from libc.math cimport fabs, pi, cos, fmin
04:
+05: import numpy as np
__pyx_t_7 = __Pyx_ImportDottedModule(__pyx_n_s_numpy, NULL); if (unlikely(!__pyx_t_7)) __PYX_ERR(0, 5, __pyx_L1_error) __Pyx_GOTREF(__pyx_t_7); if (PyDict_SetItem(__pyx_d, __pyx_n_s_np, __pyx_t_7) < 0) __PYX_ERR(0, 5, __pyx_L1_error) __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
06: cimport numpy as np
07:
08: from cython.parallel import prange
09:
+10: def apodize_edges(img: np.ndarray):
/* Python wrapper */
static PyObject *__pyx_pw_7nanopyx_4core_9transform_5edges_1apodize_edges(PyObject *__pyx_self,
#if CYTHON_METH_FASTCALL
PyObject *const *__pyx_args, Py_ssize_t __pyx_nargs, PyObject *__pyx_kwds
#else
PyObject *__pyx_args, PyObject *__pyx_kwds
#endif
); /*proto*/
static PyMethodDef __pyx_mdef_7nanopyx_4core_9transform_5edges_1apodize_edges = {"apodize_edges", (PyCFunction)(void*)(__Pyx_PyCFunction_FastCallWithKeywords)__pyx_pw_7nanopyx_4core_9transform_5edges_1apodize_edges, __Pyx_METH_FASTCALL|METH_KEYWORDS, 0};
static PyObject *__pyx_pw_7nanopyx_4core_9transform_5edges_1apodize_edges(PyObject *__pyx_self,
#if CYTHON_METH_FASTCALL
PyObject *const *__pyx_args, Py_ssize_t __pyx_nargs, PyObject *__pyx_kwds
#else
PyObject *__pyx_args, PyObject *__pyx_kwds
#endif
) {
PyArrayObject *__pyx_v_img = 0;
#if !CYTHON_METH_FASTCALL
CYTHON_UNUSED Py_ssize_t __pyx_nargs;
#endif
CYTHON_UNUSED PyObject *const *__pyx_kwvalues;
PyObject *__pyx_r = 0;
__Pyx_RefNannyDeclarations
__Pyx_RefNannySetupContext("apodize_edges (wrapper)", 0);
#if !CYTHON_METH_FASTCALL
#if CYTHON_ASSUME_SAFE_MACROS
__pyx_nargs = PyTuple_GET_SIZE(__pyx_args);
#else
__pyx_nargs = PyTuple_Size(__pyx_args);
if (unlikely((__pyx_nargs < 0))) __PYX_ERR(0, 10, __pyx_L3_error)
#endif
#endif
__pyx_kwvalues = __Pyx_KwValues_FASTCALL(__pyx_args, __pyx_nargs);
{
PyObject **__pyx_pyargnames[] = {&__pyx_n_s_img,0};
PyObject* values[1] = {0};
if (__pyx_kwds) {
Py_ssize_t kw_args;
switch (__pyx_nargs) {
case 1: values[0] = __Pyx_Arg_FASTCALL(__pyx_args, 0);
CYTHON_FALLTHROUGH;
case 0: break;
default: goto __pyx_L5_argtuple_error;
}
kw_args = __Pyx_NumKwargs_FASTCALL(__pyx_kwds);
switch (__pyx_nargs) {
case 0:
if (likely((values[0] = __Pyx_GetKwValue_FASTCALL(__pyx_kwds, __pyx_kwvalues, __pyx_n_s_img)) != 0)) {
(void)__Pyx_Arg_NewRef_FASTCALL(values[0]);
kw_args--;
}
else if (unlikely(PyErr_Occurred())) __PYX_ERR(0, 10, __pyx_L3_error)
else goto __pyx_L5_argtuple_error;
}
if (unlikely(kw_args > 0)) {
const Py_ssize_t kwd_pos_args = __pyx_nargs;
if (unlikely(__Pyx_ParseOptionalKeywords(__pyx_kwds, __pyx_kwvalues, __pyx_pyargnames, 0, values + 0, kwd_pos_args, "apodize_edges") < 0)) __PYX_ERR(0, 10, __pyx_L3_error)
}
} else if (unlikely(__pyx_nargs != 1)) {
goto __pyx_L5_argtuple_error;
} else {
values[0] = __Pyx_Arg_FASTCALL(__pyx_args, 0);
}
__pyx_v_img = ((PyArrayObject *)values[0]);
}
goto __pyx_L4_argument_unpacking_done;
__pyx_L5_argtuple_error:;
__Pyx_RaiseArgtupleInvalid("apodize_edges", 1, 1, 1, __pyx_nargs); __PYX_ERR(0, 10, __pyx_L3_error)
goto __pyx_L3_error;
__pyx_L3_error:;
{
Py_ssize_t __pyx_temp;
for (__pyx_temp=0; __pyx_temp < (Py_ssize_t)(sizeof(values)/sizeof(values[0])); ++__pyx_temp) {
__Pyx_Arg_XDECREF_FASTCALL(values[__pyx_temp]);
}
}
__Pyx_AddTraceback("nanopyx.core.transform.edges.apodize_edges", __pyx_clineno, __pyx_lineno, __pyx_filename);
__Pyx_RefNannyFinishContext();
return NULL;
__pyx_L4_argument_unpacking_done:;
if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_img), __pyx_ptype_5numpy_ndarray, 0, "img", 0))) __PYX_ERR(0, 10, __pyx_L1_error)
__pyx_r = __pyx_pf_7nanopyx_4core_9transform_5edges_apodize_edges(__pyx_self, __pyx_v_img);
int __pyx_lineno = 0;
const char *__pyx_filename = NULL;
int __pyx_clineno = 0;
/* function exit code */
goto __pyx_L0;
__pyx_L1_error:;
__pyx_r = NULL;
__pyx_L0:;
{
Py_ssize_t __pyx_temp;
for (__pyx_temp=0; __pyx_temp < (Py_ssize_t)(sizeof(values)/sizeof(values[0])); ++__pyx_temp) {
__Pyx_Arg_XDECREF_FASTCALL(values[__pyx_temp]);
}
}
__Pyx_RefNannyFinishContext();
return __pyx_r;
}
static PyObject *__pyx_pf_7nanopyx_4core_9transform_5edges_apodize_edges(CYTHON_UNUSED PyObject *__pyx_self, PyArrayObject *__pyx_v_img) {
PyObject *__pyx_r = NULL;
__Pyx_TraceDeclarations
__Pyx_TraceFrameInit(__pyx_codeobj__11)
__Pyx_RefNannyDeclarations
__Pyx_RefNannySetupContext("apodize_edges", 0);
__Pyx_TraceCall("apodize_edges", __pyx_f[0], 10, 0, __PYX_ERR(0, 10, __pyx_L1_error));
/* … */
/* function exit code */
__pyx_L1_error:;
__PYX_XCLEAR_MEMVIEW(&__pyx_t_1, 1);
__PYX_XCLEAR_MEMVIEW(&__pyx_t_2, 1);
__Pyx_XDECREF(__pyx_t_3);
__Pyx_AddTraceback("nanopyx.core.transform.edges.apodize_edges", __pyx_clineno, __pyx_lineno, __pyx_filename);
__pyx_r = NULL;
__pyx_L0:;
__Pyx_XGIVEREF(__pyx_r);
__Pyx_TraceReturn(__pyx_r, 0);
__Pyx_RefNannyFinishContext();
return __pyx_r;
}
/* … */
__pyx_tuple__23 = PyTuple_Pack(1, __pyx_n_s_img); if (unlikely(!__pyx_tuple__23)) __PYX_ERR(0, 10, __pyx_L1_error)
__Pyx_GOTREF(__pyx_tuple__23);
__Pyx_GIVEREF(__pyx_tuple__23);
/* … */
__pyx_t_7 = __Pyx_PyDict_NewPresized(1); if (unlikely(!__pyx_t_7)) __PYX_ERR(0, 10, __pyx_L1_error)
__Pyx_GOTREF(__pyx_t_7);
if (PyDict_SetItem(__pyx_t_7, __pyx_n_s_img, __pyx_kp_s_np_ndarray) < 0) __PYX_ERR(0, 10, __pyx_L1_error)
__pyx_t_4 = __Pyx_CyFunction_New(&__pyx_mdef_7nanopyx_4core_9transform_5edges_1apodize_edges, 0, __pyx_n_s_apodize_edges, NULL, __pyx_n_s_nanopyx_core_transform_edges, __pyx_d, ((PyObject *)__pyx_codeobj__11)); if (unlikely(!__pyx_t_4)) __PYX_ERR(0, 10, __pyx_L1_error)
__Pyx_GOTREF(__pyx_t_4);
__Pyx_CyFunction_SetAnnotationsDict(__pyx_t_4, __pyx_t_7);
__Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
if (PyDict_SetItem(__pyx_d, __pyx_n_s_apodize_edges, __pyx_t_4) < 0) __PYX_ERR(0, 10, __pyx_L1_error)
__Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
+11: return _apodize_edges(img)
__Pyx_XDECREF(__pyx_r); __pyx_t_1 = __Pyx_PyObject_to_MemoryviewSlice_dsds_float(((PyObject *)__pyx_v_img), PyBUF_WRITABLE); if (unlikely(!__pyx_t_1.memview)) __PYX_ERR(0, 11, __pyx_L1_error) __pyx_t_2 = __pyx_f_7nanopyx_4core_9transform_5edges__apodize_edges(__pyx_t_1); if (unlikely(!__pyx_t_2.memview)) __PYX_ERR(0, 11, __pyx_L1_error) __PYX_XCLEAR_MEMVIEW(&__pyx_t_1, 1); __pyx_t_1.memview = NULL; __pyx_t_1.data = NULL; __pyx_t_3 = __pyx_memoryview_fromslice(__pyx_t_2, 2, (PyObject *(*)(char *)) __pyx_memview_get_float, (int (*)(char *, PyObject *)) __pyx_memview_set_float, 0);; if (unlikely(!__pyx_t_3)) __PYX_ERR(0, 11, __pyx_L1_error) __Pyx_GOTREF(__pyx_t_3); __PYX_XCLEAR_MEMVIEW(&__pyx_t_2, 1); __pyx_t_2.memview = NULL; __pyx_t_2.data = NULL; __pyx_r = __pyx_t_3; __pyx_t_3 = 0; goto __pyx_L0;
12:
+13: cdef float[:, :] _apodize_edges(float[:, :] img) nogil:
static __Pyx_memviewslice __pyx_f_7nanopyx_4core_9transform_5edges__apodize_edges(__Pyx_memviewslice __pyx_v_img) {
__Pyx_memviewslice __pyx_v_pin = { 0, 0, { 0 }, { 0 }, { 0 } };
float __pyx_v_d;
float __pyx_v_c;
int __pyx_v_y_i;
int __pyx_v_x_i;
int __pyx_v_height;
int __pyx_v_width;
int __pyx_v_offset;
float __pyx_v_radius;
int __pyx_v_count;
float __pyx_v_edge_mean;
float __pyx_v_dist;
float __pyx_v_x0;
float __pyx_v_y0;
__Pyx_memviewslice __pyx_r = { 0, 0, { 0 }, { 0 }, { 0 } };
__Pyx_TraceDeclarations
__Pyx_RefNannyDeclarations
__Pyx_RefNannySetupContext("_apodize_edges", 1);
__Pyx_TraceCall("_apodize_edges", __pyx_f[0], 13, 1, __PYX_ERR(0, 13, __pyx_L1_error));
/* … */
/* function exit code */
__pyx_L1_error:;
#ifdef WITH_THREAD
__pyx_gilstate_save = __Pyx_PyGILState_Ensure();
#endif
__Pyx_XDECREF(__pyx_t_1);
__Pyx_XDECREF(__pyx_t_2);
__Pyx_XDECREF(__pyx_t_3);
__Pyx_XDECREF(__pyx_t_4);
__Pyx_XDECREF(__pyx_t_5);
__Pyx_XDECREF(__pyx_t_6);
__PYX_XCLEAR_MEMVIEW(&__pyx_t_8, 1);
__PYX_XCLEAR_MEMVIEW(&__pyx_t_18, 1);
__pyx_r.data = NULL;
__pyx_r.memview = NULL;
__Pyx_AddTraceback("nanopyx.core.transform.edges._apodize_edges", __pyx_clineno, __pyx_lineno, __pyx_filename);
goto __pyx_L2;
__pyx_L0:;
if (unlikely(!__pyx_r.memview)) {
#ifdef WITH_THREAD
PyGILState_STATE __pyx_gilstate_save = __Pyx_PyGILState_Ensure();
#endif
PyErr_SetString(PyExc_TypeError, "Memoryview return value is not initialized");
#ifdef WITH_THREAD
__Pyx_PyGILState_Release(__pyx_gilstate_save);
#endif
}
#ifdef WITH_THREAD
__Pyx_PyGILState_Release(__pyx_gilstate_save);
#endif
__pyx_L2:;
#ifdef WITH_THREAD
__pyx_gilstate_save = __Pyx_PyGILState_Ensure();
#endif
__PYX_XCLEAR_MEMVIEW(&__pyx_v_pin, 1);
__Pyx_TraceReturn(Py_None, 0);
__Pyx_RefNannyFinishContext();
#ifdef WITH_THREAD
__Pyx_PyGILState_Release(__pyx_gilstate_save);
#endif
return __pyx_r;
}
14:
15: cdef float[:] pin
16: cdef float d, c
17: cdef int y_i, x_i
+18: cdef int height = img.shape[0]
__pyx_v_height = (__pyx_v_img.shape[0]);
+19: cdef int width = img.shape[1]
__pyx_v_width = (__pyx_v_img.shape[1]);
+20: cdef int offset = 20 # number of pixels used for smooth apodization
__pyx_v_offset = 20;
+21: cdef float radius = width / 2 - offset
__pyx_v_radius = ((((long)__pyx_v_width) / 2) - __pyx_v_offset);
+22: cdef int count = 0
__pyx_v_count = 0;
+23: cdef float edge_mean = 0
__pyx_v_edge_mean = 0.0;
+24: cdef float dist = 0
__pyx_v_dist = 0.0;
+25: cdef float x0 = 0
__pyx_v_x0 = 0.0;
+26: cdef float y0 = 0
__pyx_v_y0 = 0.0;
27:
+28: with gil:
{
#ifdef WITH_THREAD
PyGILState_STATE __pyx_gilstate_save = __Pyx_PyGILState_Ensure();
#endif
/*try:*/ {
/* … */
/*finally:*/ {
/*normal exit:*/{
#ifdef WITH_THREAD
__Pyx_PyGILState_Release(__pyx_gilstate_save);
#endif
goto __pyx_L5;
}
__pyx_L4_error: {
#ifdef WITH_THREAD
__Pyx_PyGILState_Release(__pyx_gilstate_save);
#endif
goto __pyx_L1_error;
}
__pyx_L5:;
}
}
+29: pin = np.ravel(np.copy(img))
__Pyx_GetModuleGlobalName(__pyx_t_2, __pyx_n_s_np); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 29, __pyx_L4_error) __Pyx_GOTREF(__pyx_t_2); __pyx_t_3 = __Pyx_PyObject_GetAttrStr(__pyx_t_2, __pyx_n_s_ravel); if (unlikely(!__pyx_t_3)) __PYX_ERR(0, 29, __pyx_L4_error) __Pyx_GOTREF(__pyx_t_3); __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; __Pyx_GetModuleGlobalName(__pyx_t_4, __pyx_n_s_np); if (unlikely(!__pyx_t_4)) __PYX_ERR(0, 29, __pyx_L4_error) __Pyx_GOTREF(__pyx_t_4); __pyx_t_5 = __Pyx_PyObject_GetAttrStr(__pyx_t_4, __pyx_n_s_copy); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 29, __pyx_L4_error) __Pyx_GOTREF(__pyx_t_5); __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0; __pyx_t_4 = __pyx_memoryview_fromslice(__pyx_v_img, 2, (PyObject *(*)(char *)) __pyx_memview_get_float, (int (*)(char *, PyObject *)) __pyx_memview_set_float, 0);; if (unlikely(!__pyx_t_4)) __PYX_ERR(0, 29, __pyx_L4_error) __Pyx_GOTREF(__pyx_t_4); __pyx_t_6 = NULL; __pyx_t_7 = 0; #if CYTHON_UNPACK_METHODS if (unlikely(PyMethod_Check(__pyx_t_5))) { __pyx_t_6 = PyMethod_GET_SELF(__pyx_t_5); if (likely(__pyx_t_6)) { PyObject* function = PyMethod_GET_FUNCTION(__pyx_t_5); __Pyx_INCREF(__pyx_t_6); __Pyx_INCREF(function); __Pyx_DECREF_SET(__pyx_t_5, function); __pyx_t_7 = 1; } } #endif { PyObject *__pyx_callargs[2] = {__pyx_t_6, __pyx_t_4}; __pyx_t_2 = __Pyx_PyObject_FastCall(__pyx_t_5, __pyx_callargs+1-__pyx_t_7, 1+__pyx_t_7); __Pyx_XDECREF(__pyx_t_6); __pyx_t_6 = 0; __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0; if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 29, __pyx_L4_error) __Pyx_GOTREF(__pyx_t_2); __Pyx_DECREF(__pyx_t_5); __pyx_t_5 = 0; } __pyx_t_5 = NULL; __pyx_t_7 = 0; #if CYTHON_UNPACK_METHODS if (unlikely(PyMethod_Check(__pyx_t_3))) { __pyx_t_5 = PyMethod_GET_SELF(__pyx_t_3); if (likely(__pyx_t_5)) { PyObject* function = PyMethod_GET_FUNCTION(__pyx_t_3); __Pyx_INCREF(__pyx_t_5); __Pyx_INCREF(function); __Pyx_DECREF_SET(__pyx_t_3, function); __pyx_t_7 = 1; } } #endif { PyObject *__pyx_callargs[2] = {__pyx_t_5, __pyx_t_2}; __pyx_t_1 = __Pyx_PyObject_FastCall(__pyx_t_3, __pyx_callargs+1-__pyx_t_7, 1+__pyx_t_7); __Pyx_XDECREF(__pyx_t_5); __pyx_t_5 = 0; __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 29, __pyx_L4_error) __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; } __pyx_t_8 = __Pyx_PyObject_to_MemoryviewSlice_ds_float(__pyx_t_1, PyBUF_WRITABLE); if (unlikely(!__pyx_t_8.memview)) __PYX_ERR(0, 29, __pyx_L4_error) __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_v_pin = __pyx_t_8; __pyx_t_8.memview = NULL; __pyx_t_8.data = NULL; }
30:
+31: for x_i in prange(width):
{
#ifdef WITH_THREAD
PyThreadState *_save;
_save = NULL;
if (PyGILState_Check()) {
Py_UNBLOCK_THREADS
}
__Pyx_FastGIL_Remember();
#endif
/*try:*/ {
__pyx_t_7 = __pyx_v_width;
{
#if ((defined(__APPLE__) || defined(__OSX__)) && (defined(__GNUC__) && (__GNUC__ > 2 || (__GNUC__ == 2 && (__GNUC_MINOR__ > 95)))))
#undef likely
#undef unlikely
#define likely(x) (x)
#define unlikely(x) (x)
#endif
__pyx_t_10 = (__pyx_t_7 - 0 + 1 - 1/abs(1)) / 1;
if (__pyx_t_10 > 0)
{
#ifdef _OPENMP
#pragma omp parallel
#endif /* _OPENMP */
{
#ifdef _OPENMP
#pragma omp for reduction(+:__pyx_v_count) lastprivate(__pyx_v_dist) reduction(+:__pyx_v_edge_mean) firstprivate(__pyx_v_x_i) lastprivate(__pyx_v_x_i) lastprivate(__pyx_v_y_i)
#endif /* _OPENMP */
for (__pyx_t_9 = 0; __pyx_t_9 < __pyx_t_10; __pyx_t_9++){
{
__pyx_v_x_i = (int)(0 + 1 * __pyx_t_9);
/* Initialize private variables to invalid values */
__pyx_v_dist = ((float)__PYX_NAN());
__pyx_v_y_i = ((int)0xbad0bad0);
/* … */
/*finally:*/ {
/*normal exit:*/{
#ifdef WITH_THREAD
__Pyx_FastGIL_Forget();
if (_save) {
Py_BLOCK_THREADS
}
#endif
goto __pyx_L8;
}
__pyx_L8:;
}
}
+32: for y_i in range(height):
__pyx_t_11 = __pyx_v_height;
__pyx_t_12 = __pyx_t_11;
for (__pyx_t_13 = 0; __pyx_t_13 < __pyx_t_12; __pyx_t_13+=1) {
__pyx_v_y_i = __pyx_t_13;
+33: dist = (y_i - height/2)**2 + (x_i - width/2)**2
__pyx_v_dist = (__Pyx_pow_long((__pyx_v_y_i - (((long)__pyx_v_height) / 2)), 2) + __Pyx_pow_long((__pyx_v_x_i - (((long)__pyx_v_width) / 2)), 2));
+34: if dist > radius**2:
__pyx_t_14 = (__pyx_v_dist > powf(__pyx_v_radius, 2.0));
if (__pyx_t_14) {
/* … */
}
}
}
}
}
}
}
#if ((defined(__APPLE__) || defined(__OSX__)) && (defined(__GNUC__) && (__GNUC__ > 2 || (__GNUC__ == 2 && (__GNUC_MINOR__ > 95)))))
#undef likely
#undef unlikely
#define likely(x) __builtin_expect(!!(x), 1)
#define unlikely(x) __builtin_expect(!!(x), 0)
#endif
}
+35: edge_mean += pin[y_i*width + x_i]
__pyx_t_15 = ((__pyx_v_y_i * __pyx_v_width) + __pyx_v_x_i);
__pyx_v_edge_mean = (__pyx_v_edge_mean + (*((float *) ( /* dim=0 */ (__pyx_v_pin.data + __pyx_t_15 * __pyx_v_pin.strides[0]) ))));
+36: count += 1
__pyx_v_count = (__pyx_v_count + 1);
37:
+38: with gil:
{
#ifdef WITH_THREAD
PyGILState_STATE __pyx_gilstate_save = __Pyx_PyGILState_Ensure();
#endif
/*try:*/ {
/* … */
/*finally:*/ {
/*normal exit:*/{
#ifdef WITH_THREAD
__Pyx_PyGILState_Release(__pyx_gilstate_save);
#endif
goto __pyx_L20;
}
__pyx_L20:;
}
}
+39: edge_mean /= count
__pyx_v_edge_mean = (__pyx_v_edge_mean / __pyx_v_count);
}
40:
+41: for x_i in prange(width):
{
#ifdef WITH_THREAD
PyThreadState *_save;
_save = NULL;
if (PyGILState_Check()) {
Py_UNBLOCK_THREADS
}
__Pyx_FastGIL_Remember();
#endif
/*try:*/ {
__pyx_t_10 = __pyx_v_width;
{
#if ((defined(__APPLE__) || defined(__OSX__)) && (defined(__GNUC__) && (__GNUC__ > 2 || (__GNUC__ == 2 && (__GNUC_MINOR__ > 95)))))
#undef likely
#undef unlikely
#define likely(x) (x)
#define unlikely(x) (x)
#endif
__pyx_t_7 = (__pyx_t_10 - 0 + 1 - 1/abs(1)) / 1;
if (__pyx_t_7 > 0)
{
#ifdef _OPENMP
#pragma omp parallel
#endif /* _OPENMP */
{
#ifdef _OPENMP
#pragma omp for lastprivate(__pyx_v_c) lastprivate(__pyx_v_d) lastprivate(__pyx_v_x0) firstprivate(__pyx_v_x_i) lastprivate(__pyx_v_x_i) lastprivate(__pyx_v_y0) lastprivate(__pyx_v_y_i)
#endif /* _OPENMP */
for (__pyx_t_9 = 0; __pyx_t_9 < __pyx_t_7; __pyx_t_9++){
{
__pyx_v_x_i = (int)(0 + 1 * __pyx_t_9);
/* Initialize private variables to invalid values */
__pyx_v_c = ((float)__PYX_NAN());
__pyx_v_d = ((float)__PYX_NAN());
__pyx_v_x0 = ((float)__PYX_NAN());
__pyx_v_y0 = ((float)__PYX_NAN());
__pyx_v_y_i = ((int)0xbad0bad0);
/* … */
/*finally:*/ {
/*normal exit:*/{
#ifdef WITH_THREAD
__Pyx_FastGIL_Forget();
if (_save) {
Py_BLOCK_THREADS
}
#endif
goto __pyx_L23;
}
__pyx_L23:;
}
}
+42: for y_i in range(height):
__pyx_t_11 = __pyx_v_height;
__pyx_t_12 = __pyx_t_11;
for (__pyx_t_13 = 0; __pyx_t_13 < __pyx_t_12; __pyx_t_13+=1) {
__pyx_v_y_i = __pyx_t_13;
+43: x0 = fabs(x_i - width/2)
__pyx_v_x0 = fabs((__pyx_v_x_i - (((long)__pyx_v_width) / 2)));
+44: y0 = fabs(y_i - height/2)
__pyx_v_y0 = fabs((__pyx_v_y_i - (((long)__pyx_v_height) / 2)));
+45: if fabs(x0 - width/2) <= offset or fabs(y0 - height/2) <= offset:
__pyx_t_16 = (fabs((__pyx_v_x0 - (((long)__pyx_v_width) / 2))) <= __pyx_v_offset);
if (!__pyx_t_16) {
} else {
__pyx_t_14 = __pyx_t_16;
goto __pyx_L31_bool_binop_done;
}
__pyx_t_16 = (fabs((__pyx_v_y0 - (((long)__pyx_v_height) / 2))) <= __pyx_v_offset);
__pyx_t_14 = __pyx_t_16;
__pyx_L31_bool_binop_done:;
if (__pyx_t_14) {
/* … */
goto __pyx_L30;
}
+46: d = fmin(fabs(x0 - width/2), fabs(y0 - height/2))
__pyx_v_d = fmin(fabs((__pyx_v_x0 - (((long)__pyx_v_width) / 2))), fabs((__pyx_v_y0 - (((long)__pyx_v_height) / 2))));
+47: c = (cos(d * pi / offset - pi) + 1) / 2
__pyx_v_c = ((cos((((__pyx_v_d * M_PI) / ((double)__pyx_v_offset)) - M_PI)) + 1.0) / 2.0);
+48: pin[y_i*width + x_i] = <float> (c*(pin[y_i*width + x_i]-edge_mean)) + edge_mean
__pyx_t_15 = ((__pyx_v_y_i * __pyx_v_width) + __pyx_v_x_i);
__pyx_t_17 = ((__pyx_v_y_i * __pyx_v_width) + __pyx_v_x_i);
*((float *) ( /* dim=0 */ (__pyx_v_pin.data + __pyx_t_17 * __pyx_v_pin.strides[0]) )) = (((float)(__pyx_v_c * ((*((float *) ( /* dim=0 */ (__pyx_v_pin.data + __pyx_t_15 * __pyx_v_pin.strides[0]) ))) - __pyx_v_edge_mean))) + __pyx_v_edge_mean);
+49: elif fabs(x_i - width/2) > width/2 and fabs(y_i - height/2) > height/2:
__pyx_t_16 = (fabs((__pyx_v_x_i - (((long)__pyx_v_width) / 2))) > (((long)__pyx_v_width) / 2));
if (__pyx_t_16) {
} else {
__pyx_t_14 = __pyx_t_16;
goto __pyx_L33_bool_binop_done;
}
__pyx_t_16 = (fabs((__pyx_v_y_i - (((long)__pyx_v_height) / 2))) > (((long)__pyx_v_height) / 2));
__pyx_t_14 = __pyx_t_16;
__pyx_L33_bool_binop_done:;
if (__pyx_t_14) {
/* … */
}
__pyx_L30:;
}
}
}
}
}
}
#if ((defined(__APPLE__) || defined(__OSX__)) && (defined(__GNUC__) && (__GNUC__ > 2 || (__GNUC__ == 2 && (__GNUC_MINOR__ > 95)))))
#undef likely
#undef unlikely
#define likely(x) __builtin_expect(!!(x), 1)
#define unlikely(x) __builtin_expect(!!(x), 0)
#endif
}
+50: pin[y_i*width + x_i] = <float> edge_mean
__pyx_t_15 = ((__pyx_v_y_i * __pyx_v_width) + __pyx_v_x_i);
*((float *) ( /* dim=0 */ (__pyx_v_pin.data + __pyx_t_15 * __pyx_v_pin.strides[0]) )) = ((float)__pyx_v_edge_mean);
+51: with gil:
{
#ifdef WITH_THREAD
PyGILState_STATE __pyx_gilstate_save = __Pyx_PyGILState_Ensure();
#endif
/*try:*/ {
/* … */
/*finally:*/ {
__pyx_L37_return: {
#ifdef WITH_THREAD
__Pyx_PyGILState_Release(__pyx_gilstate_save);
#endif
goto __pyx_L0;
}
__pyx_L38_error: {
#ifdef WITH_THREAD
__Pyx_PyGILState_Release(__pyx_gilstate_save);
#endif
goto __pyx_L1_error;
}
}
}
+52: return np.reshape(pin, (height, width))
__Pyx_GetModuleGlobalName(__pyx_t_3, __pyx_n_s_np); if (unlikely(!__pyx_t_3)) __PYX_ERR(0, 52, __pyx_L38_error) __Pyx_GOTREF(__pyx_t_3); __pyx_t_2 = __Pyx_PyObject_GetAttrStr(__pyx_t_3, __pyx_n_s_reshape); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 52, __pyx_L38_error) __Pyx_GOTREF(__pyx_t_2); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __pyx_t_3 = __pyx_memoryview_fromslice(__pyx_v_pin, 1, (PyObject *(*)(char *)) __pyx_memview_get_float, (int (*)(char *, PyObject *)) __pyx_memview_set_float, 0);; if (unlikely(!__pyx_t_3)) __PYX_ERR(0, 52, __pyx_L38_error) __Pyx_GOTREF(__pyx_t_3); __pyx_t_5 = __Pyx_PyInt_From_int(__pyx_v_height); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 52, __pyx_L38_error) __Pyx_GOTREF(__pyx_t_5); __pyx_t_4 = __Pyx_PyInt_From_int(__pyx_v_width); if (unlikely(!__pyx_t_4)) __PYX_ERR(0, 52, __pyx_L38_error) __Pyx_GOTREF(__pyx_t_4); __pyx_t_6 = PyTuple_New(2); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 52, __pyx_L38_error) __Pyx_GOTREF(__pyx_t_6); __Pyx_GIVEREF(__pyx_t_5); if (__Pyx_PyTuple_SET_ITEM(__pyx_t_6, 0, __pyx_t_5)) __PYX_ERR(0, 52, __pyx_L38_error); __Pyx_GIVEREF(__pyx_t_4); if (__Pyx_PyTuple_SET_ITEM(__pyx_t_6, 1, __pyx_t_4)) __PYX_ERR(0, 52, __pyx_L38_error); __pyx_t_5 = 0; __pyx_t_4 = 0; __pyx_t_4 = NULL; __pyx_t_7 = 0; #if CYTHON_UNPACK_METHODS if (unlikely(PyMethod_Check(__pyx_t_2))) { __pyx_t_4 = PyMethod_GET_SELF(__pyx_t_2); if (likely(__pyx_t_4)) { PyObject* function = PyMethod_GET_FUNCTION(__pyx_t_2); __Pyx_INCREF(__pyx_t_4); __Pyx_INCREF(function); __Pyx_DECREF_SET(__pyx_t_2, function); __pyx_t_7 = 1; } } #endif { PyObject *__pyx_callargs[3] = {__pyx_t_4, __pyx_t_3, __pyx_t_6}; __pyx_t_1 = __Pyx_PyObject_FastCall(__pyx_t_2, __pyx_callargs+1-__pyx_t_7, 2+__pyx_t_7); __Pyx_XDECREF(__pyx_t_4); __pyx_t_4 = 0; __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0; if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 52, __pyx_L38_error) __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; } __pyx_t_18 = __Pyx_PyObject_to_MemoryviewSlice_dsds_float(__pyx_t_1, PyBUF_WRITABLE); if (unlikely(!__pyx_t_18.memview)) __PYX_ERR(0, 52, __pyx_L38_error) __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_r = __pyx_t_18; __pyx_t_18.memview = NULL; __pyx_t_18.data = NULL; goto __pyx_L37_return; }